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Abstract 

In this paper we study the metastability of the contact process on a random regular 
graph. We show that the extinction time of the contact process, when initialized so 
that all vertices are infected at time 0, grows exponentially with the vertex number. 
Moreover, we show that the extinction time divided by its mean converges to a unit 
exponential distribution in law. 


1 Introduction 

The contact process (C)t>o with infection parameter A on a connected, locally hnite 
graph G = {Vg, £g) is a continnous-time Markov chain that evolves as follows. For each 
t the random variable C takes valne in {snbsets of Vg}; we regard the elements of as 
infected vertices. Each infected vertex recovers with rate 1; each healthy vertex (a vertex in 
W\^t) becomes infected at rate A times the nnmber of infected neighbors. 

The contact process on a finite graph G will eventnally reach the absorbing state 0. Of 
natural interest is the time to extinction tg (dehned to be the time of the hrst visit to 0) 
when the contact process is started from the full-occupancy state. Since large connected 
graphs G will contain long linear chains, when the infection parameter A exceeds the crit¬ 
ical value Ac for the contact process on Z the contact process on G can be expected to 
survive for a long time, eventually reaching a quasi-st at ionary state that persists until, by 
chance, a large number of infected vertices almost simultaneously become healthy, leading to 
subsequent extinction. This phenomenon has become known generically as “metastability”. 
Metastability for contact process was first introduced in [3] for the one-dimensional finite 
cube. Subsequently, there have been various studies of metastability on other class of graphs, 
for example, linear chains in [18], the d—dimensional finite cube in [13], power law random 
graphs in [15] , hnite trees in [1] , and a family of hnite graphs in [IT] . 

In [7] the contact process on random d-regular graphs (d > 3) is studied. A random 
regular graph G ~ Q{n,d) is a graph chosen uniformly from all d-regular graphs with n 
vertices. Such a graph locally looks like a tree, but globally it dihers signihcantly from a 
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finite tree. Nevertheless, certain techniques and results that have been developed for studying 
the contact processes on trees will be of use here. In P HTJ HHI it is shown that when d> 3, 
there exist constants 0 < Ai(T“') < A 2 (T“') < oo that demarcate different phases for the 
contact process on an inhnite d—regular tree If A < Ai(T‘^), the contact process started 
from a hnite initial conhguration will almost surely die out; when A > A 2 (T°') it has positive 
probability of local survival; and when A is in between then it with positive probability 
survives globally but almost surely dies out locally. 

In this paper we investigate the metastable behavior of the contact process on random 
regular graphs G ~ Q{n,d) with infection parameter A > Ai(T'^). Denote by the 

contact process on G with initial conhguration A C Vg- When A = Vg we use (^t)t>o as 
shorthand; ii A = {u}, we write instead of Since the underlying graph G ^ Q{n, d) is 
random, we say that a property holds for asymptotically almost every G if the set of graphs 
in Q{n, d) which satisfy the property has probability tending to 1 as n —)■ oo. 

Throughout this paper we hx A > Ai(T'^), and we let pa > 0 be the chance that 
a contact process with initial state {0} on T'^, survives forever. To emphasize conditional 
probability and expectation given the graph G we use notations Pg and Eg- The word 
“typical” in this paper means asymptotically almost every. Also, all o(l) terms tend to 0 
uniformly in n (independent of G). 

The main results of this paper are Theorem 11.11 and II.2[ 

Theorem 1.1. There exists /3 > 0 such that for asymptotically almost every G Q{n,d), 

PG{Cxp(/3n) 7^ 0} = 1 - o(l), 

where is the contact process with initial configuration Vg- 

Theorem 1.2. For asymptotically almost every G ~ Q{n,d), the distribution o/ tg/Etg 
converges to an exponential distribution with mean 1, where tq = inf{t > 0 : = 0} is the 

extinction time of the contact process started from Vg- 

These theorems are proved in sections 2 and 3, respectively. While preparing this paper, 
we learned that J.-C. Mourrat and D. Valesin [16] have independently established Theo¬ 
rem o Because our proof is somewhat different from theirs, and because Theorem 11.11 is a 
key complement to Theorem 11.21 we include it in section 2. 

2 Exponential extinction time 

The goal of this section is to prove Theorem 11.11 

Here are some important facts. Let (fi be a contact process on with initial state {O}. 
In [m |T2| it is shown that 

Theorem 2.1. There exist constants cx,G{d) G M such that 

< E|Cf I < G(d)e""h 

Moreover, if X > Ai(T^) then cx > 0. 
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Throughout this paper ca > 0 will be the constant in the above theorem. 

In [7] (Proposition 5.2 and 5.3) it is shown that 

Proposition 2.2. Fix 0 < £ < 1/8. For asymptotically almost every G ~ Q{n,d), there are 
at least (1 — o(l))n vertices (call them “good” vertices) in G, such that for each good vertex 
u, 

^G'{'^(l+£)logn/cA ~ ~ o( 1))PA- 

Proposition 2.3. Fix 5 > 0 and 0 < £ < 1/8. Then for asymptotically almost every 
G Q{n,d), if u is a good vertex in Proposition \2.2[ then 


Pg{(1 - 5)npx < |^p+,)logn/cAl < (1 + ^)npx \ i(l+e)\ogn/c^ ^ 0} = 1 “ o(l). 


Recall that for a graph G = (Vg, V^;), the edge expansion parameter is defined as 


^E{G,k) 


min 

SCVG,|S|<fc 


\E{S,S% 

\S\ 


where E{S, S’^) C Ve is the set of edges with one vertex in S and the other vertex in It 
is shown in [5] (Theorem 4.16) that 

Theorem 2.4. Let d > 3. Then for every 5 > 0 there exists £ > 0 such that for asymptoti¬ 
cally almost every G Q{n,d), "^e^G, en) > d — 2 — 6. 


Fix an integer M > 0. Suppose ?7 C Vg is of size an, where a > 0. We remove every 
vertex in U whose M-neighborhood is not a tree, and denote the remaining vertex set by U'. 
We claim = an — o{n). Here we are using the following fact shown in [10] (Lemma 3.2): 

Proposition 2.5. For asymptotically almost every G ~ Q{n,d), it has at most o{n) vertices 
whose \ \og^_i\ogn\-neighhorhoods in G are not tree-like. 

Since the cardinality of U' and U are on the same order of magnitude, without loss of 
generality, let us assume that all vertices in U have tree-like M-neighborhoods in G. 

We classify vertices in U into 2 categories by looking at their M-neighborhoods in G 
in the following way. For v E U, let B{v,M) be the induced subgraph containing all 
vertices in n’s M-neighborhood in G. B{v, M)\{v} has d connected components, call them 
Gi{v), G 2 {v), , Gd{v). If Gi{v) contains no other vertices in U, call it a free branch of depth 

M of V. 


• Color V black if at least one of Gi{v), G 2 {v),..., Gd{v) is a free branch of depth M of 

V. 


• Color V white if none of Gi{v), G 2 {v),..., Gd{v) is a free branch of depth M of n . 

Proposition 2.6. Fix M G N. There exists e = s{M) > 0, such that for asymptotically 
almost every G ~ Q{n,d) the following statement holds: for any set U G Vg satisfying 
\U\ < en and that every vertex in U has its M-neighborhood in G being a tree, then U has 
at least |t/|/4 black vertices. 
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Proof. We will construct a subset of vertices W C Vq- First of all, W contains all vertices 
in U. Moreover, we are going to add some vertices into W based on the white vertices of 
U. Let V E U he a. white vertex. In each of Ci{v), C' 2 (n),..., Cd{v) there must be at least 
another vertex in U. Suppose x E U D Ci{v), then for the pair {v,x), we add into W every 
vertex along the (unique) geodesic between v and x. We repeat this operation for every 
possible pair {v,x) to obtain W. 

Such constructed W contains 3 types of vertices: black vertices of U, white vertices of 
U, and the vertices which are added by the above procedure (color them grey). Now let us 
count their contributions to E{W, W^). 

• A white vertex will contribute 0 edge to E{W, W^). This is because all of its d neigh¬ 
boring vertices are already in W by our construction. 

• A black vertex can contribute at most d edges to E(W, W^), possibly fewer. 

• A grey vertex can contribute at most d — 2 edges to E(W, possibly fewer. This 
is because by our construction a grey vertex must be sitting on the geodesic between 
two other vertices in U and therefore at least 2 out of its d neighboring vertices are 
already in W. 


Suppose in U there are w white vertices, b black vertices. Then g, the number of grey vertices 
in W, satishes g < {d + d{d — 1) -1- ■ ■ ■ -|- d{d — := Nmw. 

Due to Theorem 12.61 there exists Sm such that on a typical random regular graph G, 
^EiG^EMn) > d — 2 — {3d — 8)/{3Nm + 4). This forces the following inequality (provided 
b + w + g< SMn), 

OJ _ O 

Ow + db + {d — 2)g > E(W., W^) > {d — 2 -—- 7 )('^ + b + g). 

3Nm + 4 


Together with g < Nmw, we conclude that 

b 1 

- > -• 

b + w 4 


Therefore take = EmUNm + 1) (this guarantees that if b + w < e'j^^n then b + w + g < 
Em^) . As long as \U\ < E'j^n and every vertex in U has its M-neighborhood in G being a 
tree, then U has at least |t/|/4 black vertices. □ 


We need the following result about the growth rate of the severed contact process on a 
tree, shown in [7] (Proposition 2.2). A severed contact process {?7P}t>o is a version of contact 
process on with initial conhguration {0} where we do not allow infections to come across 
d — 1 edges connected to O. 

Proposition 2.7. There exists A = A(A, d) > 0 such that > Aexp(cAt), for all t > 0. 

Let Am he the following hnite graph. Fix the root O in T'^. We hrst remove d — 1 
edges connected to O in T'^. For the remaining graph, the connected component of the 
M-neighborhood of O is called A^. In A^, O and vertices at distance M from O are of 
degree 1; other vertices are of degree d. Am is a tree. 

Let be a contact process on Am with initial conhguration {O}. The following 

corollary is easily obtained from Proposition 12.71 
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Corollary 2.8. For every N > 0, there exist T > 0 and M G N such that W\ri^'^^\ > N. 
Now we can prove Theorem II.1[ 

Proof of Theorem li.il Let G be a typical graph as in Proposition 12.4[ 12.51 and 12.61 

In Corollary 12.81 take iV = 10 so that there exist T and M satisfying > 10. 

Withont loss of generality assnme T > 1. Fnrthermore, take L > 0 large enongh snch that 
Emin(|p^’^*^|, L) > 9. Fix snch choices of T, M and L. 

Let e = e{2M) be the constant in Proposition 12.61 As long as f/ C Vg is of size en, 
then there will be at least en/2 vertices in U snch that each vertex has its 2M-neighborhood 
being a tree, and by Proposition 12.61 there will be at least en/% black vertices. We ennmerate 
the black vertices to be ni, ^2 • • • Wfc where k > en/S. 

Now for each n*, it has one free branch of depth 2M (if more than one, specify one). Add 
Vi (and the edge connected to Vi) to its free branch of depth M (as a snbgraph of the branch 
of depth 2M) and we obtain a snbgraph isomorphic to Am- So each Vi is associated with a 
copy of Am- Here we nse Am instead of A 2 M to ensnre they are disjoint. For each Vj, we 
rnn an independent contact process on its copy of Am, call it {rii*}t>o- In a standard way 
we can conple and together so that is always dominated by . 

Let Xi = min(|r 7 ^*|, L), then EX* > 9, and 0 < X* < L. Fnrthermore (Xj)i<j<fc are i.i.d. 
random variables, so by Hoeffding’s ineqnality, 

P{^(A', - 9) < -k} < exp (-g) < exp [-^) . 

i=l k / 

in other words, after time T, with probability at least 1 — exp(—£n/(4L^)), we will observe 
at least 8/c > en infections in 

To snmmarize, as long as we rnn the contact process with initial conhgnration whose 
cardinality is en, then after time T, with probability more than 1 — exp(—£n/(4L^)), the 
ontcome will have cardinality at least en. Let (3 = £/(8L^), we conclnde 

^ 0} > 1 - exp(/3n) exp j = 1 - o(l)- 


□ 

Remark 2.9. One can slightly change the above proof to show the following statement. 
There exists a constant eo > 0. For every 0 < e < eo, there exists /3^ > 0, such that for 
asymptotically almost every G ~ Q{n,d), for any U gVq with \U\ = en, we have 

PG{5e'ip(&„) 0} = 1 - 0 ( 1 ). 


3 Metastability 

Let tg = inf{t > 0 : = 0} be the extinction time of the contact process started from fnll 

occnpancy on G. In this section we are devoted to proving Theorem 11.21 The key ingredient 
is the following proposition in [13] (Proposition 2.1): 
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Proposition 3.1. Suppose there is a sequence of graphs Gn = (V„,£^„). Let Tn be the 
extinction of the contact process started from full occupancy on Gn- Also, for arbitrary 
U C Vn, we couple fff and in the standard way. Suppose there exist two sequences of 
positive real numbers a{n) and h{n), both tending to infinity, such that as n ^ oo, 

1 . a{n)/b{n) 0 ; 

S- ^ 0 } ^ 1 - 

Then r„/Er„ converges in distribution to an exponential distribution with mean 1. 

Now from Theorem 11.11 we can take b{n) = exp(/3n). It suffices to find an appropriate 
sequence a{n) such that item 1 and 2 in Proposition 13.11 hold. 

Lemma 3.2. There exists 7 > 0, such that for asymptotically almost every G ^ Q{n,d), for 
any two vertices u, v of G, 

PG{t’ e 52“ > n-r 

Proof. By [2], the diameter of a typical random regular graph is (l + o(l)) log^_^ n. Therefore 
on such a graph, for any pair of vertices {u,v), the graph distance between u and v is no 
more than (1 + o(l)) log^_;^ n. Assume dist{u,v) = I, this means we can hnd a sequence of 
vertices (wi)o<i<z, with wq = u, wi = v, and that Wi is connected to tCj+i in the graph. 

One way of observing {n G ^ n) follows. In time interval [f, i+1] for 0 < i <1—1, 

we require the infection at vertex Wi to go to tCj+i, and stay there alive till the end of the time 
interval. This will happen with probability p > {1 — e~^)e~^. In time interval [/, 2 log^_i n], 
we require the infection at v to stay alive. This will have probability at least 
Therefore, the overall probability is at least p’-e~^‘^^°^d-in-i) > (1 — ng- 2 iogrf_in > 

n~'^ where 7 = 2 /\og{d — 1 ) + 21og(e^/(e'^ — l))/log((i — 1 ). □ 

Lemma 3.3. There exists 60 > 0, such that for asymptotically almost every G ~ Q{n,d), 
for any v eVg, 

Pg{IC 2 t| > <5on I ^ 0 } = 1 - o(n"^). 

Proof. For asymptotically almost every G ~ Q{n,d), there are n — o{n) good vertices 
by Proposition 12.21 In particular, there is at least one good vertex. We £x one good 
vertex in G, call it w. Our idea is as follows. We chop the time interval [0,n^'>'] into 
[0,T], [T,2T],..., [(M - 1)T,MT], where M = [n^yT\ and T to be specified. We will 
construct an event and see if it happens in each interval [iT, {i + 1)T]. The chance that it 
happens in each time interval is on the order of n~'^. This event is constructed such that as 
long as in at least one of these intervals the event happens, then with probability approaching 
1 at the end we will see S^n infections. 

Here are more details. Let T = 2 log^_i n + (1 + e) logn/c\, where this £ > 0 is the same 
as in Proposition 12.21 and 12.31 We say that we observe a success in time interval [(f — 1)T, iT] 
if the following two events happen. 

1. W E ^p_i)T+21ogrf_in ; 
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2 . 


> npx/2. 


Conditional on Q 27 7 ^ 0, for sure 7 ^ 0, so by Lemma [Sill ^ ^(i-i)r+ 2 iogrf_i n) 

pens with probability at least n~'^ (notice that the event {w G 'Cp_i)r+ 2 iog^ ^nl positively 
correlated with the event {^^(27 7 ^ 0})- Given {w G ^(i_i)T+ 2 iogrf jnl; since ta is a good vertex, 
by Proposition 12.31 ^ 'npx/2 will happen with probability at least px/2. 

Therefore overall we have order (n^'^/logn) trials, each with success probability at least 
Pxn~'^/2, so it is easy to conclude that the chance of having at least 1 success is 1 — o(n“^). 
Given that we observe a success, which means that we observe pxn/2 infections at some time 
between [0,71^"^], from the proof of Theorem 11.11 and Remark 12.91 we know that the chance 
of pxn/2 infections not lasting exponentially long time before the size of infections shrinks 
to 5on is exponentially small in n, where (5o can be taken as min(£o,PA/2). Therefore the 
overall probability is (1 — o(n“^)) x (1 — o(e“^’^)) = 1 — o(n“^). □ 

We call a sequence of vertices vo,...,vl G Vg a path of length L if UjUj+i G £g for 
0 < z < L — 1. We say such path has endpoints vq and vl- A path of length 0 is allowed, 
where the endpoints are identical. 

Lemma 3.4. For every r > 0, for asymptotically almost every G ~ Q{n,d), there exist 
Lr eN and Cr > 0, such that for any U,W C Vg with \U\ = \W\ = rn, there exits CrU paths 
in G so that 


• Each path is of length no more than Lj.. 


• Each path has one endpoint in U and the other endpoint in W. 

• None of the paths share the same vertex. 

Proof. From [1] and [6], there exits > 0 so that for a typical d-random regular graph G its 
(vertex) isoperimetric constant is at least hd. Without loss of generality, we assume r < 1/2. 
For z G N, Let t/, = {n G Vg : dist{v, U) < 1 } and hfo = {n G Vg : dist{v, W) < i}. By the 
definition of the isoperimetric constant, Ui has cardinality at least (1 + hd)rn, and that Ui+i 
has cardinality at least (1 + hd)\Ui\, provided that \Ui\ <n/2. The same holds for hfo. 

Take /{2r))'\ +2, we claim that |foL^/ 2 | > 2 d+ 2 ferf ^' because from 

the way we choose for sure there will be some k < Lr/2 — 1 so that k is the smallest integer 
satisfying \Uk\ > n/2. If \Uk\ > 2 d+ 2 hd ^ done. Otherwise, since < zz/ 2 , the 

number of vertices in Uk that are adjacent to Uf is at least hd\Uf\, and that the number of 
vertices in Uf that are adjacent to Uk is at least hd\U^\/d. In this case 

It/wil > ir^i + L|c7j|, 

Since \Uk\ > n /2 and |17^| > n — f^i^n, we have 


1 ^. , n hd 

t/fc+i > - + — 
' -2d 



d + 2hd \ 
2d + 2hd ) 


d T 2hd 

—;-n. 

2d + 2hd 
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Now we are able to find > 2 d+ 2 hd ^i ''^^ere k,j < Lr/2. By the inclusion- 

exclusion principle, 

\Uknw,\> 


-n. 


d-\- hd 

Now each vertex n G f/^ fl Wj naturally corresponds to a path of length no more than 
Lr/2 -|- Lr/2 = Lj. with one endpoint in U and the other in W. Notice that there exits 
Ml^ G N, so that any path of length no more than Lr can intersect at most other paths 
of length no more than in G, therefore we can pick at least 


L- 


hd 


-n\ 


' Ml^ d-\- hd 

non-intersecting paths which satisfy all requirements stated in the lemma. 


□ 


Proposition 3.5. Let a{n) = 2n^'^ -|- Ls^, where ho is as in Lemma [Ql and Lsq is as in 
Lemma fS^ Then for asymptotically almost every G Q{n,d), 


sup ^ 0,C(n) 7^ Qn)} = o(l). 

UcVg 


Proof. It suffices to show 


sup Pg{C,„| # 0.C(„) # CT„)} = »(’> ). 

veVc 

because for arbitrary U gVg, 

PelCtn) 0.«h) 7^ d5)} < El’GfQ") 7^ 7^ «?«)>■ 

vGU 


( 1 ) 


Moreover, in order to show ([T]) it suffices to show the following inequality by a union 
bound, 

sup PciQn) ^ 0, w ^ C(n),u G = o(n"2). (2) 

ve.Va,uGVa 

Recall the graphical representation and the dual contact process in [9]. To show ([2]), we 
chop the time interval [0,a(n)] into [0,71^"^ -|- Ls^] and -|- Lso,a{n)]. We run for time 

+ Lsq in the first subinterval and run the dual process Q for time in the second 
sub interval, where time t for the dual process corresponds to time a{n) — t for the original 
process. In particular, time 0 for the dual process corresponds to time a{n) for the original 
one, and time for the dual process corresponds to time + Ls^ for the original one. 

Notice that and are two independent processes. Then we have the following upper 

bound of ([2]), 


PcIQn) 7^ 07 « ^ Qn)7 W e 7^ 07 Ot 7^ 07 C^-r+Lso ^ = 0} 

<^ g { C 2 i + Ls ^ n ^“27 = 0 I 7^ 0, ^ 2 ^ 7^ 0}. 

Our goal is to show the following bound holds uniformly for u, v: 

n o, = 01 o, ^ 0, o, # 0} = o(n-^). (3) 









By Lemma 13.31 

Pg{|&,I > 6on, 10,1 > Son \ O, 7^ 0,0, # 0} = (1 - oirr^yf = 1 - o(n-^). (4) 

Given {|0,l > 0on, I0,l > ^ 0 ^}) by Lemma [3.41 there will be cs^n non-intersecting paths 
of length no more than Lsg with one endpoint in ^^27 and the other in ^“ 27 - We call it a 
success on a path vo,Vi,... ,Vk with endpoints Vq G .^^27 and Vk G (^“ 27 , if in time Ls^, the 
infection at Vq spreads to within this path. For non-intersecting paths, successes on paths 
are independent events. Also, using similar argument as in proof of Lemma [3.31 it is easy to 
see that there exits pso > 0 so that the probability of a success is at least ps^ on a path of 
length no more than 

In order to observe {^^ 27 _,_i^ H ^“27 7 ^ 0}, it suffices that on one of these cs^n paths we 
observe a success, each having probability at least ps^ >0. Therefore the chance of observing 
at least 1 success is at least 1 — o(n“^) by an easy binomial calculation. 

Combining all arguments above, we see 

^g{C^^+Lsq bl C 27 7^ 0 I Czt 7^ 0, In27 7^ 0} 

>Pg{| 0 .| > 5on, IC 27 I > ^on | O. ^ 0,67 ^ 0 } x (1 - o{n-^)) 

=l-o(n- 2 ), 

which hnishes our proof. □ 

Proof of Theorem I i. M By Proposition 13.11 take b{n) = exp(/3n) as in Theorem 11.11 and 
a{n) = 2'n?'^ -|- as in Proposition 13.51 □ 
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